Nonlinear guided waves and symmetry breaking in left-handed waveguides 



m 
o 
o 

» 

o 
O 



' 

o 1 

> ■ 

Oh- 

q ; 
o ■ 

• i-H . 

C0 . 

>V 

43 ' 



> 

(N 
O 

o 
o 

: 

o . 

o : 

>\ 
43 : 

Oh. 



13 



Ilya V. Shadrivov 

Nonlinear Physics Group, Research School of Physical Sciences and Engineering, 
Australian National University, Canberra ACT 0200, Australia 

We analyze nonlinear guided waves in a planar waveguide made of a left-handed material sur- 
rounded by a Kerr-like nonlinear dielectric. We predict that such a waveguide can support fast and 
slow symmetric and antisymmetric nonlinear modes. We study the symmetry breaking bifurcation 
and asymmetric modes in such a symmetric structure. We analyze nonlinear dispersion properties 
of the guided waves, and show that the modes can be both forward and backward. 
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Experimental realization of composite structures with 
simultaneously negative dielectric permittivity and neg- 
ative magnetic permeability \\ revived the theoretical 
interest to the left-handed materials (LHMs). Such ma- 
terials were first considered theoretically by Veselago , 
who predicted negative refraction at the interface be- 
tween LHM and conventional dielectric, inverse Doppler 
and inverse Vavilov-Cherenkov effects. Although LHMs 
do not exist in nature, artificial metallic composites p] 
can effectively be described as materials having both di- 
electric permittivity and magnetic permeability negative 
in the microwave frequency range. Nowadays, the efforts 
have been made to design LHMs for optical frequencies 
using metallic nanowires 0. Additionally, it was shown 
that photonic crystals demonstrate negative refraction 
under some conditions 0, , thus resembling the most 
famous feature of LHMs. 

In I960's Veselago predicted that a slab of LHM with 
some particular parameters surrounded by usual dielec- 
tric can focus light emitted by a point source located at 
one side of the slab to the point image on the other side of 
the slab Q Later, Pendry suggested @ that such a lens 
can focus not only propagating waves, but also evanes- 
cent waves, thus creating a perfect image of the source. 
Although, hot debates on the possibility to create such 
a perfect lens followed this suggestion (see, e.g. 0), it 
is clear that the planar structures containing LHMs have 
a high potential for the future image processing appli- 
cations. In particular, it was already shown, that the 
periodic structure of alternating layers of negative refrac- 
tion material and dielectric possesses an unusual band- 
gap structure 0, |j| , when the average refractive index on 
the structural period is equal to zero, and it can be used 
for the complex beam shaping . 

Recent studies of the waveguiding by the structures 
with LHMs have shown that the properties of guided 
modes in such systems differ essentially from those of 
conventional waveguides. So far, linear waves guided by 
an interface and by a slab of LHM (see, e.g. as 
well as nonlinear surface waves |12| have been studied. It 
has been shown that the mode structure and frequency 
dispersion may have an unconventional form and unex- 
pected properties. 

The study of linear waves of a LH slab waveguide 



has shown some of their peculiar properties, such as the 
absence of the fundamental mode, mode double degen- 
eracy, and the existence of both forward and backward 
waves. In this paper, we study nonlinear guided modes in 
a waveguide formed by a slab of linear LHM embedded 
into nonlinear dielectric, and we show that symmetric, 
antisymmetric and asymmetric waves are supported by 
such a waveguide, and we study their properties. We 
predict, that in nonlinear regime additional modes which 
do not exist in linear problem appear. We demonstrate 
that the propagation of the wavefronts (characterized by 
the phase velocity) with respect to the direction of the 
energy flow (Poynting vector) depends on the propaga- 
tion constant, and the waves can be both forward and 
backward travelling. 

To study nonlinear guided waves in a nonlinear waveg- 
uide, we consider a LH slab with real negative dielectric 
permittivity €2 and real negative magnetic permeability 
\ii surrounded by a nonlinear dielectric [see inset in Fig. 
n^a)] with constant magnetic permeability fix and dielec- 
tric permittivity e\ depending on the intensity of the elec- 
tric field 



ei {\E\ 



€1 + a\E\ 



(1) 



where we assume a > which corresponds to the 
sclf-focusing nonlinear medium. The waveguide is uni- 
form along y-axis. To be specific, we consider TE- 
polarized guided waves, which are governed by the non- 
linear Helmholtz equation 



AE y + k 2 e(x,\E y \ 2 )fi(x)E y = 0, 



(2) 



where A = d 2 /dx 2 + d 2 /dz 2 is a two-dimensional Lapla- 
cian, e (x, \E y \ 2 ^ and \i (x) are the dielectric permittiv- 
ity and magnetic permeability, ko = uj/c is a free-space 
wavenumber, 10 is an angular frequency, and c is the speed 
of light. The guided waves in a similar system, but with 
a conventional dielectric core have been studied analyti- 
cally in Ref. [13. 

We look for stationary guided modes in the form 
E = ty(x) exp (ihz), where the transverse mode struc- 
ture ^S(x) can be determined from the equation 



d 2 V 

dx 2 



[k 2 e (x, |*| 2 ) - h 2 ] * = 0. 



(3) 
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FIG. 1: (color online) (a) Dependence of the normalized 
power of guided modes p on the propagation constant 7. Pa- 
rameters are: L = 2, ei = 1, yLti = 1, €2 = —2, /i2 = —2. 
Vertical line 7 = 2 separates the fast (to the left of the line) 
and slow modes. Bold solid curve - symmetric mode, dotted 
- antisymmetric, dashed - asymmetric, thin solid line - power 
of two solitons in nonlinear medium vs. propagation constant. 
Dashed rectangle is magnified below in (b). Numbers indicate 
parameters for which the mode structure is shown in Fig. [5] 



We introduce the dimensionless variables (x',z') = 
(koX,koz), 7 = h/ko, ip — tyy/fiia. For notational sim- 
plicity we omit the primes below, and in the dimension- 
less variables Eq. © can be rewritten in the form 



+ [c 2 /i2 - 7 2 ] i> = 0, \x\ <L, 
+ [ei/xi - 7 2 ] V> + M V = 0, \x\>L, (4) 



dSp_ 
dx 2 



where 2L is a dimensionless thickness of the slab. 

At this point, we separate the guided waves into fast 
and slow modes. The fast modes have the phase velocity 
larger than the phase velocity of light in a homogeneous 
medium of the core. For such modes, 7 2 < £2/^2, and 
their localization is caused by the total internal reflec- 
tion of light from the cladding, resembling localization 
of waves in a dielectric waveguide. For the slow modes, 




FIG. 2: (color online) Mode structure calculated for the pa- 
rameters indicated in Figs. Q^a,b). 



7 2 < £2/^2, and the wave guiding resembles localization 
of the surface waves. 

Solutions of Eq. JIJ for guided modes can be found in 
the form 



V / 2Kisech[Ki(x — xi)], x < —L, 
ip(x) = { Asm(k 2 x) + B cos(k 2 x), \x\ < L, 
V / 2Kisech[Ki(x — X2)], x > L, 



(5) 



where k 2 = 7 2 — e\H\, k\ — €211-2 — 7 2 , and A, B, x\,X2 
are constants determined from continuity of the tangen- 
tial components of the electric and magnetic fields at the 
interfaces at x = —L and x = L. The fast modes cor- 
respond to k\ > 0, while for the slow modes fcf < 0. 
Solutions (JSJ) at |jc| > L have the form of sech- functions 
or solitons which are centered at x — x\ and x — X2 at 
either side of the slab. The modes with x\ < —L and 
X2 > L have the field maxima at the corresponding side 
of the waveguide. From the linear theory, it follows that 
solutions for the stationary modes can be found sepa- 
rately for the symmetric and antisymmetric modes. As 
has already been shown |13j . even in a symmetric non- 
linear dielectric waveguide asymmetric modes can exist. 
To find the asymmetric waves general solution (JSJ of the 
differential equation Q should be considered here. Equa- 
tions for the parameters x\, X2 of the guided modes can 
be found from the boundary conditions by eliminating 
the constants A and B, 



sech(Mi) [/3tanh (Mi) - tan (k 2 L)] 



sech(M 2 ) [tan(fc 2 L) 
sech(Mi) [/3tanh(Mi) 



■/3tanh(M 2 )] 
- tan" 1 (k 2 L)] 



sech(M 2 ) [tan -1 (k 2 L) + /3tanh(M 2 )] 



1, 



= 1, 



(6) 



where Mi = K\{L + x\), M 2 = K\{L — x 2 ), and (3 = 
KijU 2 /fc 2 /ii. Corresponding equations for the symmetric 
and antisymmetric modes (for which —x\ = .x 2 = xq) 
can be obtained in the form: 



(3 tanh [n\{L — xq)] — ± tan ±:L (fc 2 L), 



(7) 
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FIG. 3: (color online) (a) Dependence of the normalized 
power of guided modes p on the propagation constant 7. Pa- 
rameters are: L = 2, ei = 1, fj,i = 3, £2 = —2, /X2 = —2. 
Vertical line 7 = 2 separates the fast (to the left of the line) 
and slow modes. Bold solid curve - symmetric mode, dotted 
- antisymmetric, dashed - asymmetric, two thin solid lines - 
dispersion of one (lower curve) and two (upper curve) soliton 
states in the corresponding nonlinear media. Numbers indi- 
cate the parameters for which the mode structure is shown in 
(b). 



where (+) corresponds to the symmetric modes, while 
(— ) to the antisymmetric ones. Solution of Eqs. JJJJ can 
be found numerically after reducing them to a single 
equation, or the analytical solution can be obtained in 
the way discussed in Ref. 01, where the problem has 
been solved separately in two different cases, AB > 
and AB < and the explicit expressions for unknown 
parameters have been obtained. 

Energy flow in a stationary guided mode has the only 
component along the waveguide and it can be found as 
an integral of the Poynting vector: 



P 



c 

47T 



E y H x dx. 



(8) 



In the normalized form, the energy flow can be written 
as 

p = — P = -^1 [2- tanh (Mi) - tanh (M 2 )l + 

c 2 fix 



7 

2fi 2 



L (B 2 + A 2 ) + Sin {2k2L) (B 2 - A 2 ^ 



2k 2 



(9) 



Note, that the energy flow in the mode can be both pos- 
itive and negative with respect to the wavevector, since 
inside the slab the energy propagates in the opposite di- 
rection to that outside the slab, and the total energy flow 
depends on a ratio of these quantities. 

Dependence of the normalized energy flow p on the di- 
mensionless wave number 7 for both fast and slow modes 
is shown in Fig. ^ Parameters in the figure are chosen 
in such a way that only one fast mode with a single node 
(close to the point 4) exists in the linear case (at low in- 
tensities; for the solution of linear problem and for the 
choice of parameters see, e.g. Ref. [ill). Thin solid curve 
shows the power of two solitons in a homogeneous nonlin- 
ear cladding versus the propagation number. The mode 
structures are shown in Fig. [21 where each plot demon- 
strates the transverse wave profile corresponding to the 
numbered point in Fig. ^ 

The modes with the parameters close to the dotted 
curve (points 1, 5, 8 in Fig. ^) resembles two in-phase 
and out-of-phase solitons at either side of the waveguide. 
Closer to the dotted curve the soliton centers in the non- 
linear media move further away from the waveguide. The 
symmetry breaking bifurcation appears on the symmet- 
ric mode branch. The asymmetric mode dispersion curve 
in the fast region ends in the symmetric mode branch [see 
Fig. ^ (a)] , while in the slow wave region the asymmet- 
ric mode disappears, when the amplitude of the wave at 
one interface becomes the same as the amplitude of the 
soliton on the other side of the waveguide. Two struc- 
tures (points 6,7 in Fig. [5J obtained at the same value 
of the propagation constant show the point of the sym- 
metry breaking, where a slight asymmetry can be seen 
in the mode structure shown in the example 7 in Fig. [21 
Symmetric and asymmetric modes can be both forward 
and backward travelling (p > and p < 0, respectively). 

For the parameters indicated in Fig. ^ only one fast 
guided mode with a single node exists in the linear limit, 
while in the nonlinear regime the modes with zero, one 
and two nodes appear (see Fig. [21 . With increasing the 
waveguide thickness, one more symmetry breaking point 
appears on the antisymmetric mode branch, when the 
slab parameter L exceeds some threshold value. More- 
over, more high-order modes can be supported by the 
structure, and the nonlinear dispersion diagram becomes 
more complicated. 

Nonlinear dispersion shown in Fig. [HIa) is obtained for 
the parameters when two slow modes exist in the linear 
case (at low intensities) . The transverse structure of the 
modes, corresponding to the numbered points in Fig.[3fa) 
is shown in Fig.[3fb). For symmetric and antisymmetric 
modes dispersion curves start from zero energy points 
corresponding to linear modes in the slow-wave region 
of the plot, and they end at the curve representing the 
power of two solitons when the modes have the structure 
of two in-phase or out-of-phase solitons shifted from the 
waveguide at either side of it. Nonlinearity can only re- 
duce the propagation constant of these modes (bold solid 
and dotted lines), thus increasing the wave phase veloc- 
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ity. These symmetric and antisymmetric modes become 
fast with the growth of the intensity. This becomes possi- 
ble due to the nonlinear change of the index of refraction 
profile in the waveguide cladding. This index modifica- 
tion is caused by a high field amplitude in the part of the 
cladding next to the core, as one can see from the plots 
(1,3,4,6) in Fig. El 

The symmetry breaking bifurcation point is located 
at the symmetric mode branch. The asymmetric mode 
branch ends at the curve corresponding to the power of 
one-soliton, when the mode is represented by a soliton 
infinitely shifted from the waveguide from one side of it 
and a vanishingly small field trapped inside the waveg- 
uide (point 9 in Fig. EJ). Plots (4,5) in Figs. E{b) shows 
the symmetric and asymmetric modes structure for the 
same value of the propagation constant at the point of 
the symmetry breaking. 

The two presented dispersion characteristics of the 
nonlinear waves in a LH waveguide surrounded by a Kerr- 
like nonlinear medium show the general dispersion prop- 



erties of the low-order fast and slow modes, which are 
qualitatively similar in other parameter regions. 

In conclusion, we have studied nonlinear guided waves 
in a waveguide structure created by a slab of linear LH 
material surrounded by a nonlinear conventional dielec- 
tric. We have shown that fast and slow symmetric, asym- 
metric and antisymmetric modes can exist in this sys- 
tem. On the dispersion diagram, the symmetry break- 
ing occurs for both symmetric and antisymmetric mode 
branches. In the latter case, there exists a threshold value 
of the slab thickness above which the asymmetric modes 
appear. The guided waves can be both forward and back- 
ward travelling. While the direction of the energy flow is 
determined by a source, the type of the wave determines 
the phase velocity direction, i.e. the phase front propa- 
gates from the source in the forward waves, and to the 
source, in the backward waves. 

The author thanks Yuri S. Kivshar and Andrey A. 
Sukhorukov for useful discussions. This work was par- 
tially supported by the Australian Research Council. 



[1] D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat- 

Nasser, and S. Schultz, Phys. Rev. Lett. 84, 4184 (2000); 

M. Bayindir, K. Aydin, E. Ozbay, P. Markos, and C. M. 

Soukoulis, Appl. Phys. Lett. 81, 120 (2002). 
[2] V. G. Veselago, Usp. Fiz. Nauk 8, 2854 (1967) [Sov. Phys. 

Usp. 10, 509 (1968)]. 
[3] V. A. Podolskiy, A. K. Sarychev, and V. M. Shalaev, 

Journ. of Nonlin. Opt. Phys. and Materials 11, 65 (2002) 
[4] M. Notomi, Phys. Rev. B. 62, 10696 (2000). 
[5] H. Kosaka, T. Kawashima, A. Tomita, M. Notomi, T. 

Tamamura, T. Sato, and S. Kawakami, Phys. Rev. B 58, 

10096 (1998). 
[6] J. B. Pendry, Phys. Rev. Lett. 85, 3966 (2000). 
[7] N. Garcia and M. Nieto-Vesperinas, Phys. Rev. Lett. 88, 

207403 (2002); J. B. Pendry, Phys. Rev. Lett. 91, 099701 

(2003); N. Garcia and M. Nieto-Vesperinas, Phys. Rev. 

Lett. 91, 099702 (2003); D. R. Smith, D. Schurig, M. 



Rosenbluth, S. Schultz, S. A. Ramakrishna, and J. B. 

Pendry, Appl. Phys. Lett. 82, 1506 (2003); 
[8] J. Li, L. Zhou, C. T. Chan, and P. Sheng, Phys. Rev. 

Lett. 90, 083901 (2003); 
[9] I. V. Shadrivov, A. A. Sukhorukov, and Yu. S. Kivshar, 

Appl. Phys. Lett. 82, 3820 (2003). 
[10] R. Ruppin, Phys. Lett. A 277, 61 (2000); T. A. Leskova, 

A. A. Maradudin, and I. Simonsen, SPIE Proc. 4447, 6 

(2001). 

[11] I. V. Shadrivov, A. A. Sukhorukov, Yu. S. Kivshar, Phys. 

Rev. E 67, 057602-4 (2003). 
[12] I. V. Shadrivov, A. A. Sukhorukov, Yu. S. Kivshar, 

A. A. Zharov, A. D. Boardman, and P. Egan, 

arXiv:physics/0305126 
[13] N. N. Akhmediev, Zh. Eksp. Teor. Fiz. 83, 545 (1982) 

[Sov. Phys. JETP 56, 299 (1982)]. 



